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Abstract 

We study the decays of r _ — > v T P~^ [P = ir~ ,K~) in the light front quark model. 
We calculate the form factors and use them to evaluate the decay widths. We find 
that, in the standard model, the decay widths are 1.62 x 10~ 2 (3.86 x 10~ 3 ) T T -_ tUrlT - 
and 1.91 x 1(T 3 (5.38 x 10~ 4 ) Y t -^ VtK - with the cuts of £ 7 = 50 (400) MeV and t = 
800 (1200) MeV for r~ — > z/ T 7r~7 and r~ — > v T K~ r f, respectively. We also show that 
with including the radiative decay widths, the experimental rate for r _ — > v T P~ can be 
explained. 



1 Introduction 



Tau is the only charged lepton which decays into hadrons. Theoretically, the hadronic 
t decays can provide us with valuable information on strong interaction. The simplest 
rare radiative r decays are t~ — > z/ T P _ 7 with P being the pseudoscalar mesons of n 
and K. The contributions to the decays can be divided into "internal-bremsstrahlung" 
(IB) and "structure-dependent" (SD) parts in terms of the photon emission. For the IB 
contribution the photon emits from the r lepton as well as external hadron, while the SD 
one from intermediate states described by the vector and axial- vector form factors, Fy t A- 
Our main task is to calculate these form factors which are functions of t 2 = (p + q) 2 where 
q (p) is the four momentum of 7 (P) and Mp < t 2 < m 2 . In general, the momentum 
dependences of F V} a could also help us to determine the bound state wave functions of 
mesons. 

In this paper, we will use the light front quark model (LFQM) to evaluate the matrix 
elements in r — > v T P^i decays. The LFQM has been widely applied to study the form 
factors of weak decays [1, 2, 3, 4, 5, 6, 7]. It is the relativistic quark model [8] in which 
a consistent and relativistic treatment of quark spins and the center-of-mass motion can 
be carried out. Moreover, the meson state of the definite spins could be relativistically 
constructed by the Melosh transformation [9]. They are many advantages in the LFQM. 
For example, the light front wave function is manifestly Lorentz invariant as it is expressed 
in terms of the momentum fraction variable (in "+" component) in analogy with the 
parton distributions in the infinite momentum frame. The kinematic subgroup of the 
light front formalism has the maximum number of interaction-free generators including 
the boost operator which describes the center-of-mass motion of the bound state [8]. 

The paper is organized as follows. In Sec. 2, we present the matrix elements and 
study the form factors in the P — > 7 transitions within the framework of the LFQM. We 
calculate the decay widths of r — > u T Pj in Sec. 3. We also compare our results with those 
in literature [10, 11]. We give our conclusions in Sec. 4. 

2 Form Factors in the Light Front Quark Model 

Similar to the radiative meson decay, the decay amplitude for 

r-(l) -> u T (k)P-(pHq) (1) 
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with P = 7r or K can be written as [10, 12, 13] 

M = M IB + M SD , (2) 



M-ib = Gf cos 9 c ef P m T u(k)(l + 75) 



111 _|_ ^ _ ill 
p • q 21 ■ q I ■ q 



v(l), 



ka G F cos0 c e f . „ Fy 

D = 7y| U^vpaL M e VP CT — 

+«(fc)(l + 7fe)[(p-?V-(c-pV]«(0^}, (3) 

where U 1 = -u(/c)7 M (l — j 5 )u(l) and FA,y are the form-factors corresponding to the vector 
and axial-vector currents, defined by 

F 

(l(q),P(p)\qi^ 5 b\0} = e-^- [(p • g)e* - (e* • p)gj , 

( 7 (g),P(p)|g 7 ^|0) = ie^, ra ^ w ? V, (4) 

with e M being the photon polarization vector and g (p) the four momentum of 7(F), 
respectively. The form factors F^y in Eq. (4) depend on t 2 = (p + q) 2 for which the 
allowed range is Mp <t 2 < m 2 , in contrast to < t 2 < Mp in P + — > l + vi r ). 

2.1 Light front formalism 

To calculate of the hadronic matrix elements in Eq. (4), one usually lets t + = to have a 
spacelike momentum transfer. However, since the momentum transfer should be always 
timelike in a real decay process, in this work we use the frame of t± = with the physically 
accessible kinematic region of < t 2 < t^ ax . Within the light front formalism, the meson 
bound state, which consists of a quark q 1 and an anti-quark q 2 with the total momentum 
p and spin S, can be written as 

\M(p,S,S z )) = Jldhjidk^nY^ip-h-h) 

x 2$ ra -(fc 1 ,fc 2 ,A 1 ,A 2 )6+(fc 1 ,A 1 )d+(fc 2 ,A 2 )|0>, (5) 

Ai A2 

where 

dk + dk 



[dk] ~ 2{2-Kfk+ ' 

{b\>(k'), b\(k)} = {dy(k'),d[(k)} = 2(2vr) 3 S 3 (k' - k) 5 yx , (6) 

and fci(2) is the on-mass shell light front momentum of the internal quark (71(172) ■ The light 
front relative momentum variables (x, k±) are defined by 

kf = xip + , k 2 = x 2 p + , X1 + X2 — I, 

h± = xip_L + kj_, k 2 ± = x 2 p± - k± . (7) 



In the momentum space, the wave function Q SSz is given by 



® SSz (ki, k 2 , Ai, A 2 ) = R^ 2 (x, k ± ) <f>(x, k ± ), 



ss z 



(8) 



where <f>(x, k±) describes the momentum distribution amplitude of the constituents in the 
bound state and Rff\ 2 constructs a spin state (S, S z ) out of light front helicity eigenstates 
(A1A2), expressed by 
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R S x %(x,k ± ) = E( A il^i/( 1 -^^' m i)l s i)( A 2|^L(a;,-A; ± ,m 2 )| S2 )(- Sl - S2 | 1 S^). (9) 



S1,S2 



In Eq. (9), are the Pauli spinors and 1Z M is the Melosh transformation operator, 
given by 



TZ M (x, k±,mi) 



mi + XjM + ia ■ k± x n 



with 



Ml 



'{mi + XiM ) 2 + k 



m i + k± m 2 + k\ 



Xl 



x 2 



n = (0, 0, 1). 



(10) 



(11) 



We note that Eq. (9) can, in fact, be expressed as a covariant form [1] 



where 



and 



R S x%(x,k ± ) 



hvi ft, 2 



y/2 M 



u(ki, Xi)rv(k 2 , A 2 ), , 



M = y jM*-{m 1 -m 2 )< 



u(k, X)u(k, A) = 

A 



m+ fit ^ x . ., , , m— , 
-j^, 'Ev(k,X)v(k 1 X) = — j ^ 

A 



(12) 



(13) 



with 



r = 75 

r = - ${s x ) -\ 



(pseudoscalar, 5 = 0), 

g ■ (fci ~ jfe) 
M + mi + m 2 



e"(0) 



—e ± (±l)-p ± ,0,e ± (±l) 
p+ 

1 ^-M 2 +pi + \ 



7; 



+ ,P ,P± 



(vector, 5 = 1), 



ei(±l) = =F(l,±i 



(14) 



(15) 
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The normalization condition of the meson state is given by 

(M(p>, S', S' z )\M(p, S, S z )} = 2(2vr) *p+P(f - p) Ws»5. , (16) 

which leads to 

/w w *'* J|2=1 - (17) 

In principle, the momentum distribution amplitude <f>(x, k±) can be obtained by solving 
the light front QCD bound state equation [14, 15]. However, before such first-principle 
solutions are available, we would have to be contented with phenomenological amplitudes. 
One example that has been often used in the literature is the Gaussian type wave function: 



^ dfv-< ( k 



«*.*>■) = -£3-1' (18) 



where N = A(ir / uj^) ^ , k = (k±, k z ), and dk z /dx = eie 2 /x(l — x)M with k z being defined 
through 



x = — — 1 — x — 62 — ^ z , ti = Jmf + k 2 (19) 

ei + e 2 ei + e 2 



by 



In particular, with appropriate parameters, the wave function in Eq. (18) describes sat- 
isfactorily the pion elastic form factor up to t 2 ~ 10 GeV 2 [16]. 



2.2 The Form Factors of F A and F } 



v 



The one-loop diagrams that contribute to Fv,a are illustrated in Figure 1. From the 
figures, the hadronic matrix elements in Eq. (4) are found to be 

(7(?), J P(p)l^(l-7fe)Q|0> = / 

i(V; + m Q ) i(-y; + m g ) 7 + m q ) 



:{ta.(1- 75)^ 



i 2 — + A; 2 2 — + 9 A; 2 — m 2 + ie 
-(q^QX^K)} (21) 

where Ap is a vertex function related to Qq bound state of the meson P, k 2 = p — k[ and 
k± — t — k[ — k 2 + q- After integrating over the LF momentum k~{ in Eq. (21), we get 

(7(<7 7 ),nP)l^(l-75)Q|0> 

= /Vfeill ,- \- (in k '-) h .-\- -(q^QX^K)}, (22) 
Jo Ui -k lon k ion k 2 -k 2on J 
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where 



[d 3 k[] dktdk 



2(2n)% + k'+kf ' 
7 ^lfe-„ = Tr [l^ 1 -75)(Vi + m Q )~f 5 {-tf 2 +m q )ie q 4(tfi +m 9 )J , 

^ion = * ) ^1(2) = Pon ~~ ^2(l)on ' ^1 = ^2on + 9 7 o" ' 

with {on} representing the on-shell particles. By considering the "good" component 
pL — +, the meson wave function in Eq. (8) and Melosh transformation in Eq. (12) are 
related to the bound state vertex function A P , given by [1, 17]: 



fcj.) . (24) 



^2 ~~ ^2on \/2 M 

Moreover, the matrix elements in Eq. (4) become 

(l(q),P(p)\dl + l 5 Q\0) = e^-(f? ± -q ± ), 

(l(q),P(p)\qi + Q\0) = -ie^-^qj. (25) 

Here we have used the LF monentum variables (x, k±) and worked in the frame that the 
transverse momentum is purely longitudinal, i.e., t± = 0. We note that t 2 = t + t~ > 
covers the entire range of momentum transfers. Thus, the relevant quark variables for 
Figure. 1 are 

k+ = (1 + x)q + , k 2 = xq + , k 1± = (1 + x)q ± - k±, k 2± = xq± - k± . 
ki~ = (1 — x')p + , k 2 = x'p + , k l± = —(1 — x')q± + k ± , k 2i _ = —x'q± — k ± , (26) 

where x' (x) is the momentum fraction of the anti-quark in the meson (photon) state. At 
the quark loop, it requires that 

^2(1) = ^2~; ^2(1)± = , (27) 

for Figure la(b). Therefore, the trace in Eq. (23) can be easily carried out. The form 

factors Fa and Fy in Eq. (25) are then found to be 

p ,.2s AM /* dxd 2 k ± j 2(l + 2x)(x + x') m q -x'(m Q -m q )k 2 ± &<S>(x',k 2 ± )^ 
= -4M Pj / o ^ff^f 8 ) 

(1 + 2x){l - x" + x) m Q + (1 - x")(m Q - m q )k\Q $ (x", k]_) \ 



3 tjiq + k\ x"{l - x") 



2 f dxd 2 k ± j 2(x + x') (1 + 2x)m q - Ak]Q $ (x', k\) 

F v (t) - -4M P J Q 2(27r) 3 e j 3 ^(i_ x >) 

(1 - x" + x) (1 + 2x)m Q - fffcje $ (x", fcj) | 

3 m^ + Ai x"(l-x")j' 1 j 



respectively, where 



p+ M 2 
q+ t 2 - Ml ' 



v 

A = x'(m,Q — m q ) — 2xm q , 
B = —(l—x")(mQ — m q ) — 2xm q , 

»r/ X(l-X) \ 1/2 /rffcT / fc 2 



e 



^2(M 2 - (m Q - m,) 2 ) / V dx \ 2uj 
1 d$(x,k 
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M 



$(x,kl) dk 2 ± ' 
/t 2 -M 2 \ /t 2 -M 2 \ 

x 1 = x(-^yx"=i-x(-^y k ={k^k z ). (so 

3 Decay widths 

To compute numerical values of the form factors, we use /„- = 0.925, itlq = m q = m u = 
ma = 0.25, M n = 0.14, u n = 0.3 for the 7r meson, and fx = 0.113, itlq = m s = 0.4, 
m q = m>u = 0.25, Mk = 0.495, ujk = 0.37 for the K meson in GeV [18, 19], respectively. 
We start with the decay of r — > u T n^. We define x — 21 • q/m 2 and y — 21 ■ p/m 2 . In the 
r rest frame, x (y) corresponds to the photon (pion ) energy of E^^, expressed in units 
of m T /2 as 

2Ey 2E W 

x = and y = . (31) 

m T m T 

In terms of x and y, one has the following kinematics: 

2 

TfL 

p-q = -^(x + y-l-r), 

t 2 = ( s - k) 2 = (p + q) 2 =m 2 T (x + y- 1) . (32) 

The physical allowed regions for x and y are given by: 

< x < 1 — r , 

\-x+-H— <y< 1+r, (33) 
1 — x 



with 

/,; V (i.21 x 10 ;! (34) 



.m 7 

where we have used m T = 1.777 GeV. We now calculate the differential decay rate of 
t~ — > z/ T 7T _ 7, given by 

dr{T ^ = _L^> (1 - q - p - . (35) 



In the tau rest frame, one has 



d 2 T 



m 



T -,\M\\ 



where 



dx dy 2567T 3 



\M\ 2 = \M IB \ 2 + \M SD \ 2 + 2Re(M IB M* SD ) 



(36) 



(37) 



By writing the decay width T in terms of the three different source as in Eq. (38), i.e., 



total 



— Tib + T,sd + T/jVT 



(38) 



and the non-radiative decay (r — > v t ti ) width 

G 2 F \V ud \ 2 c f 2 



in 



m A r (l - r) 2 = 2.44 x l(T iU MeV , 



(39) 



we obtain that 



d 2 T 



IB 



dx dy 

SD 



d 2 T 



a mi 



dx dy 



\F v \ 2 Pvv(x,y) 



16tt f 2 ml 

+2Re(F v FX)p V A(x,y) + \F A \ 2 p AA (x, y) 



v 

(1-0 



d 2 T 



INT 



a mt 



dx dy 



[Pintv(x, y) Re(F v ) + pi NT a(x, y) Re(F A )} 



where 

PiB(x,y) 

Pvv(x,y) 
Paa(x,v) 
PvA(x,y) 
piNTA(x,y) 

piNTv(x,y) 



^Um^—^ a '^ v ' • — — — (1 _ r)2 

[r 2 (x + 2) - 2r(x + y) + (x + y - l)(x 2 - 3x + 2 + xy)]{r -y + l) 

(r — x — y + l) 2 x 2 
-[r 2 (x + y) + 2r(l - y)(x + y) + (x + y - l)(x 2 + y 2 - x - y)\ , 

Pvv(x,y,r) , 

[r 2 (x + y) + (1 — x — y)(y — x)](r 2 — x — y + 1) , 
(r — x — y + 1) (r — y + 1) 



, (40) 



x 



[r 2 - 2r(x + y) + (l-x + y)(x + y - l)](r - y + 1) 
(r — rr — y + l)x 

To simplify our calculations, we now introduce A as a new parameter 



(41) 



A = —r = x + y- 1. 



mt 



(42) 
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The kinematical boundaries for x and A are given by: 



X-r <x< 1 - -, 
A 

r < A < 1 . 



(43) 



By integrating the variable x in the phase space, from Eq. (41), we derive the expressions 
of differential decay widths for the invariant mass spectrum as: 
dT IB 



dX 



a 
2^ 



(1- A)(r 2 + 2rA-4A + A 2 ) 



dT 



SD 



dX 

dTiNTV 

dX 

dTjNTA 

dX 



+ (r 2 A + 2rA - 2A - 2A 2 + A 3 ) In A 

a m\ (A-l) 2 (A-r) 3 (l + 2A) 
48^ 7>f J? 

a m 2 T (A-r) 2 (l-A + AlnA) 



A 2 -rA(l-r) 



2 ' 



,2 



a 



ml 



+ {rX-2X- A 2 ) In A 



A 

1-A)(r-2A-1) 

" A 



(\Fv 
Re(Fy) 



+ \F A 



r 

2\ T— *V T -K 



(1-r) 2 ' 



A 



(44) 



Note that for the IB part in Eqs. (40) and (44) the contribution is infrared divergent 
when x closed to and, moreover, there is also an enhancement in the limit A — > r. This 
means that the IB term contains the logarithm divergent as t 2 — > m 2 . To obtain the 
decay width of r — > z/ r 7T7, a cut on the photon energy is needed. The differential decay 
width dTis/dX as a function of A = t 2 /m 2 in terms of T T ^ UTn is shown in Figures 2 and 
3 with two different cuts of E 1 = 50 and 400 MeV, respectively. In Table 1, we list the 
integrated decay width ratio of R n = T T ^ UT7T7 /T T ^ VTn for the two cuts. As shown in Table 

Table 1: Integrated Decay width ratio for r — > v t t\^(. 



Integrated decay width ratio 


IB 


SD 


INT 


Sum 


Ref.[10] 


Ref.[ll] 


10- 3 i^(£ 7 > 50) MeV 


13.1 


1.48 


1.61 


16.2 




14.8 


10- 3 i? w (£ 7 > 400) MeV 


1.48 


1.48 


0.90 


3.86 


1.16 


2.76 



1, our results for the decay width of r — > u t -k^ are larger than those in Refs. [10] and 
[11], respectively, due to our bigger SD contribution. From Table 1, it is interesting 

to see that T^r which depends on the photon energy cut can be as large as T$d- 

For the decay of r — > v T K^^ we use T t ^ Vt k = 1.759 x 10 -11 MeV. The differential 

decays width for A = t 2 /m 2 > 0.2 GeV in terms of F t _, UtK is plotted in Figure 4. From 
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the figure, we see that the SD contribution to the decay width is dominant for values of 
t 2 . The reason is that the large kaon mass suppresses the IB contribution. In Table 2, we 
show the decay width ratio of Rk = T t ^>v t Ki/Tt^>v t k with two different invariant mass 
cuts of t = 800 and 1200 MeV. Form Table 2, we find that the results in Ref. [11] are 



Table 2: Integrated decay width ratio for r — > v T K~f. 



Integrated Decay width ratio 


IB 


SD 


INT 


Sum 


Ref. [11] 


10~ 4 R K {t > 800 MeV ) 


6.51 


7.36 


5.26 


19.13 


35.8 


10- 4 i^(t > 1200 MeV ) 


0.72 


3.80 


0.86 


5.38 


9.1 



twice as large as our predictions. 
We now study the ratio 

R = ^kBt{t^uP) 
Br[r — > ev e v T ) 

and examine both theoretical and experimental values. In the standard model, the ratio 
is given by: 

R SM = 0.646 , (46) 

while the recent experimental average is [18] 

R ex P = o. 66 ioo ± 0.00725 . (47) 

Clearly, there is a discrepancy between R SM and R exp . However, it is believed that it 
arises from the radiative corrections. At 0(a), the radiative corrected decay width is 
found to be 

r(r _ -> i/ T 7r(7)) = r(r _ -> i/ T 7r _ ) + r(r _ -> z/ r 7r- 7 ) ~ 2.48 x 10~ 10 MeV (48) 

with Ej > 50 MeV. Similarly, for the mode with K, one has that 

r(r _ -> i/ T A""(7)) ~ 1.61 x 10" 11 MeV (49) 

with the same photon energy cut as the pion mode. In the standard model, the radiative 
corrected width of r — > evv [20] is given by 

r(r _ -»• e _ i/ e P r (7)) = (4.033 ± 0.005) x 10" 10 MeV. (50) 

From the Eq. (48)- (50), we obtain 

R Theor,0(a) _ 0.655 (51) 

which agrees with experimental data in Eq. (47) within the errors. 
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4 Conclusions 



We have studied the decays of r — > v T 7r(K)~/ in the light front quark model. We 
have calculated the form factors and used them to evaluate the decay widths. We have 
found that, in the standard model, the decay widths are 1.62 x 10~ 2 (3.86 x 10 -3 ) F T ^ UTn 
and 1.91 x 1(T 3 (5.38 x 1(T 4 ) T t ^ UtK with the cuts of £ 7 = 50 (£ 7 = 400) MeV and 
t = 800 (1200) MeV for r — > u t tt^ and r — > v T K^/, respectively. We have also shown 
that with including the radiative decay widths, the experimental rate for r — > u T {j^K) 
can be understood. In future, as the tau-charm factories will produce a large number of 
samples of r lepton pairs, these rare r decays should be precisely measured and thus the 
structure-dependent form-factors can be well tested. In this case, the decays can be used 
to probe new physics similar to Ki 3 [21] and decays [22]. 
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t k' 2 7r(p) 

(b) 



Figure 1: Loop diagrams that contribute r — > u T Tvy. 



13 




Figure 2: Differential decay width of T(r — > z/ T 7T7) in terms of T t ^ Ut7T as a function of 
A = t 2 /ml with E 1 > 50 MeV. The dot, dash, dash-dot and solid curves stand for the 
contributions of IB, SD, INT, and total parts, respectively. 
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Figure 3: Same as Figure 2 but with £ 7 > 400 MeV. 
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